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Introduction 


NASA has initiated the development of quiet supersonic business 
Jets 

Return of commercial supersonic flight will allow passengers to travel 
over the continental U.S. within hours and complete international 
business trips within a single day. 

NASA awarded contract for preliminary design of a low boom flight 
demonstrator for Quiet Supersonic Technology project (QueSST) 
Most efforts of the design are focused on reducing the sonic boom 
ground signature, however the noise constraints during takeoff and 
landing at subsonic speeds must be satisfied. 

Computational aeroacoustic (CAA) tools can be used to assess the 
new designs at lower speeds. 

This work represents the first part of a systematic validation effort in 
jet noise prediction capability for NASA Ames Launch Ascend and 
Vehicle Aerodynamics Code (LAVA). 


Outline 


o Experimental Setup 


Experimental Setup 


o Small Hot Jet Acoustic Rig (SHJAR), which is located in the 
Aeroacoustics Propulsion Lab (AAPL) at NASA Glenn 
Research Center 


Perspective 


view of PIV 
SHJAR — measurement 
kins device 


Bridges et. al. (NASA-TM-2011-216807) 


Experimental Setup 


o Baseline axisymmetric convergent Small Metal Chevron 
(SMCOOO0) nozzle at Set Point 7 (SP7) 
o Nozzle axis in downstream flow direction is marked as 180° 


Bridges et. al. (NASA-TM-2011-216807) SP7 
Acoustic Mach number U,.,/c 0.9 
Jet temperature ratio T./T 0.835 
Nozzle pressure ratio NPR 1.861 
Nozzle Diameter D rane 
30° Reynold number Re, 1 Mio 
X,/D= 0.242 [-] ‘ Reynolds number Re 800 
. Boundary layer thickness 0.0128 D 
oa.) ee . similar to: Bres et. al. (AIAA-2015-2535) 
SMCO0O0 teen “Bruit et vent” jet-noise facility at 


Nozzle oe —_— Centre d’Etudes Aerodynamique et Termique 
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o Computational Methodology 


Computational Methodology 


LAVA Framework (Kiris et al. Aerospace Science and Technology, Volume 55, 2016) 
o Computational Fluid Dynamics Solvers 

¢ Cartesian, Curvilinear, and Unstructured Grid Types 

¢ Overset Grid and Immersed Boundary Methods 

e Steady and Unsteady RANS (Reynolds Averaged Navier-Stokes) 

¢ Hybrid RANS/LES (Large Eddy Simulation), LES and LBM Capabilities 
o Acoustic Solver 

° Linear Helmholtz Scattering Code 

° Permeable Surface Ffowcs Williams-Hawkings Propagation 
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Cartesian Immersed Boundary — Unstructured Arbitrary Polyhedral Overset Structured Curvilinear 


Computational Methodology 
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Computational Methodology 


3-D Structured Curvilinear Overset Grid Solver 

o Spalart-Allmaras turbulence model (baseline turbulence model) 

Low-Dissipation Finite Difference Method (Housman et al. AIAA-2016-2963) 

o 6th-order Hybrid Weighted Compact Nonlinear Scheme (HWCNS) 

o Numerical flux is a modified Roe scheme 

o 6t/5' order blended central/upwind biased left and right state 
interpolation 

o 2°¢_order accurate differencing used for time discretization 

Hybrid RANS/LES Models 

o Delayed Detached Eddy Simulation (DDES) model with modified 
length scale (Housman et al. AIAA-2017-0640) 

o Zonal RANS-NLES (numerical LES) with user selected zones of URANS, 
NLES, and wall-distance based hybrid RANS-NLES (see paper for details) 

Synthetic Eddy Method 

o Coupling Methodology between RANS and LES to introduce realistic 

turbulent eddies (arrin et al. int. Journal of Heat and Fluid Flow 30 ) 9 


Computational Methodology 


o When transitioning from RANS to LES in wall-bounded flows it is 
necessary to insert meaningful three-dimensional content at the 
interface 

o The synthetic eddy method (SEM) is one approach which adds eddies 
in such away that first and second order turbulent statistics can be 
Satisfied. (approx. from the RANS solution with Bradshaw hypothesis) 


Jet Case SP 7 
AxSEM = X exit — XsemM . = . 
54.6 < Axspy < 556 


Vi 4.010 -1,008 0.007 4.005 0.004 0.002 4.001 0.001 0.002 O04 0.005 0.007 0008 0010 


Computational Methodology 


o Unsteady RANS until jet is fully developed and 
eddy viscosity maximum has plateaued 
Ar=1-10"(s1;0-41s} Restart simulation with Hybrid RANS/LES Models 

| until transient behavior washed out 
o Ignore transients which are taken at first 30000 
time-steps and restart simulation 
RANS/LES o Record Volume data at 100 kHz sampling 


frequency for greater than 0.02 seconds 


At = 1-10~° [s] ; nt > 30000 


(approx. 205 convective time units) 
ieee cowie 

Processors 1392 (has) 260 (ivy) 960 (has) 

final Hybrid Wall-Clock Time [day] 12.5 

RANS/LES Sub-iterations 5 
At = 1-107 [s] Convergence 2-4 orders every sub-iteration 
Stra. = 16.25 , Ste, = 0.008 
Teonv © 205 Number Eddies (SEM) = 5000 5000 
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o Structured Overset Grid System 


Structured Overset Grid System 


Baseline (256 M) 

Coarse ( 28M) 

Refined (106 M) 

Seven point overlap 

No orphan points 

Minimum stencil quality 0.9 
Baseline follows Bogey et. al 
(AIAA-2016-0261) 
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Structured Overset Grid System 
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-o Circumferential refinement in 


axial and radial direction 
neal 


Bres et. al. (AIAA-2015-2535) 
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o Circumferential refinement in 
axial and radial direction 
Bres et. al. (AIAA-2015-2535) 
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Structured Overset Grid System 
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o Circumferential refinement in 
axial and radial direction 
Bres et. al. (AIAA-2015-2535) 
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-0.7 321.15 
core -0.7 0.50 
wall 0.0 34.50 
core 0.0 0.06 
shear 0.5-—25.0 10.50 
core Oey 225.40 1.05 


circumferential/radial AR 


wall -0.7 436.82 
core -0.7 1.00 
wall 0.0 221.00 
core 0.0 1.00 


shear 0.5-25.0 1134 
core 0.5-—25.0 1.00 


Structured Overset Grid System 
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dx/BL_ delta: 0.00 0.08 0.15 0.23 0.31 0.38 0.46 0.54 0.62 0.69 0.77 0.85 0.92 1.00 
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dtheta/BL_delta: 0.00.1 0.20.20.3 0.40.5 0.5 0.6 0.7 0.8 0.8 0.9 1.0 
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Computational Results 


Flow Field Visualization: lso-contour of Q-criteria colored by axial velocity 


_U/Ujet: 0.00 0.08 0.15 0.23 0.31 0.38 0.46 0.54 0.62 0.69 0.77 0.85 0.92 1.00 


How can we 


improve/remove 
the 2D structures? 


Eo en ace 
Visible 2D structures —m a EES 
near the nozzle exit. .-— 
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Computational Results 


Indicator Function: DDES-256M 


o Indicator function f, indicates if in 
RANS or LES mode. 

o Stays in RANS mode in nozzle 
interior and quickly transitions to 
LES downstream of nozzle lip 

o Retains large eddy viscosity 
throughout the boundary layer 


fd: 0.00.1 0.20.20.30.40.5 0.5 0.6 0.7 0.8 0.8 0.9 1.0 


Ur / Loo: DDES-256M 


o Shielding function RANS-NLES: 


1 
{f= 1-= 5 [1 = tanh(€y(dyau i dy))] 


dwvat 2: Walldistance 
do : transition distance (user) 
€,  : blending (user) 


Eddy Viscosi ty Ratio: 0 4 8 121519232731 3538 42 46 50 


Computational Results wasal 


Indicator Function: DDES-256M Indicator Function: RANS-NLES-SEM-1L06M | 


fd: 0.00.1 0.20.20.30.40.5 0.5 0.6 0.7 0.8 0.8 0.9 1.0 fd: 0.00.10.20.20.30.4050.50.60.70.80.80.9 1.0 


Lr / oo. DDES-256M Lr / loo. RANS-NLES-SEM-106M 
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Eddy Viscosity Ratio: 0 4 8 121519232731 3538 42 4650 Eddy Viscosity Ratio: 0 4 8 121519232731 3538 42 46 50 


Waffle Cone 
Structure 


U/Ujet: 0.00 0.08 0.15 0.23 0.31 0.38 0.46 0.54 0.62 0.69 0.77 0.85 0.92 1.00 


RANS-NLES-SEM-28M 


U/Ujet: 0.00 0.08 0.15 0.23 0.31 0.38 0.46 0.54 0.62 0.69 0.77 0.85 0.92 1.00 


DDES-256M 


U/Ujet: 0.00 0.08 0.15 0.23 0.31 0.38 0.46 0.54 0.62 0.69 0.77 0.85 0.92 1.00 U/Ujet: 0.00 0.08 0.15 0.23 0.31 0.38 0.46 0.54 0.62 0.69 0.77 0.85 0.92 1.00 2 1 


RANS-NLES-106M | RANS-NLES-SEM-106M 


Computational Results asa 


RANS- — SEM dalle Mesh 5 Quasi-2D waffle cone 


structures at nozzle exit 

Size of turbulent structures 
appears to be too large inside 
nozzle 

Structures deep in the 
boundary layer show very little 
azimuthal variation 


Features are elongated and 
too highly correlated in both 
the streamwise and azimuthal 
direction 

Do we have realistic, fully 
developed BL turbulence at 
exit? 
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U/Ujet: 0.20 0.22 0.23 0.25 0.26 0.28 0.29 0.31 0.32 0.34 0.35 0.37 0.38 0.40 


Computational Results — Near-Field@oy 


Near-Field Comparison 


Centerline 


o Near field turbulent statistics computed from DDES, RANS- 
NLES and RANS-NLES-SEM models for comparison with PIV 
data from the SHJAR 

o Comparison of measurements to data at lip-line (z/R=1) and 
Centerline (z/R=0) 
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Near-Field Comparison: Time-Averaged Centerline 


Computational Results — Near-Field 
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Near-Field Comparison: Time-Averaged Lip-line 


NAS: : 


Computational Results — Near-Field 
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Computational Results — Near-Field 
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Computational Results — Far-Field S 
FWH Permeable Surface 


o Generated surface 
triangulation imbedded 
within the overset grid 

o FWH surface spans the 
entire axial domain of the 
computational grid 

o Edge length of the 
triangles set to 5mm 


Observers 100 D away 


Computational Results — Far-Field S 
FWH Permeable Surface 


o Generated surface 
triangulation imbedded 
within the overset grid 

o FWH surface spans the 
entire axial domain of the 
computational grid 

o Edge length of the 
triangles set to 5mm 


under consideration 


Observers 100 D away 


Computational Results — Far-Field €& 


o Volume solution is interpolated to the FWH surface at a sampling rate 


of At = 0.00001 s (100 kHz) 


o Total time sample is split into 5 windows (or segments) with 50% 


O 


overlap at an St,,, = 0.02 


Integrands of the FWH permeable surface formulation are computed 
over each window independently 
© Q,,F1,F,,F; 
¢ Hanning Window is applied in the time-domain 
e FFT is applied and stored for computing far-field observer noise levels 
FWH surface integrals are computed for each observer over each 
window independently 
¢ 360 observers, uniformly distributed along the azimuth, are generated for 
each jet axis angle (60°,90°, 120°, 150°) 
e¢ The PSD is ensemble averaged over the 360 observers and the PSD is 
multiplied by sqrt(8/3) to recover the RMS levels lost from Hanning Window 


Finally, the PSD spectrum is averaged over the 5 windows for the final comparison 
to the experimental consensus. 29 
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Far-Field Comparison: PSD Spectrum at 100D from exit a0 


Computational Results — Far-Field (3 
Far-Field Comparison: Band-Limited OASPL (0.08 < St <8.0) 
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Computational Results — Far-Field Gy} 
Far-Field Comparison: Band-Limited OASPL (0.08 < St < 8.0) 
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Computational Results — Far-Field 


Time-Domain Pressure Associated with Peak Frequency (1100Hz) in 150° 
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Summary 


The hybrid RANS/LES approach, within the LAVA framework, using 
structured curvilinear overlapping grids has been applied to the 
prediction of jet noise and compared to existing near-field PIV and 
far-field microphone data. 

Demonstrated improvements: 


¢ Hybrid RANS-NLES reduces the delay in transition to 3D turbulent structures 
and improved lip-line RMS prediction 

¢ SEM eliminates delay even further 

Completed far-field acoustic propagation 

¢ Mach wave radiation noise in the jet direction is well-captured 

e Sideline noise caused by turbulent fluctuations is over-predicted, likely do to 
elevated lip-line RMS at nozzle exit 

BL needs to be resolved better inside of nozzle for further 

improvements 


Future Work 


a Posterior! Error Analysis - AIAA-2017-0978 Anisotropic grid-adaptation 
Axial (j-dir) E in LES of wall-bounded and free shear flows, 


| Toosi and Larsson 
—o Analyze the difference in turbulent 


kinetic energy using the resolved 
velocity field with filtered version of 
resolved velocity field. 

o Independent filtering in each direction 
leads to anisotropic measure for 


resJ: 0 2 4 6 8 101213151719 21 2325 


= refinement 


Radial (I-dir) = Circumferential (k-dir) 
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Future Work 


a Posteriori Error auc ees AIAA-2017-0978 Anisotropic grid-adaptation 
Axial (j-dir) ~~ __ in LES of wall-bounded and free shear flows, 


Toosi and Larsson 


o Resolution in streamwise direction is 
lacking the most. 
o The error estimate has largest 
"magnitude in circumferential direction. 
SSS _o Radial direction pretty good. 
——_—_—— ——— -o Improved mesh (191 M) for further 
|, a |S __investigation of SP7 and all SP3 runs 


Circumferential (k-dir) Circumferential (k-dir) 
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Future Work asa 
LES with explicit subgrid-scale (SGS) model and SEM 


o No RANS downstream of SEM 
location 

o Waffle cone structures inside 
nozzle reduced 

o Artificial turbulence from SEM 
decays towards nozzle exit due 
to lack of resolution 

o Recommended resolution: 
wall-resolved As7._.. = 20 
(12.5k points) 
wall-modeled As;,;,- = 0.16 
(2450 points) 


We need wall-modeled LES in 
order to resolve the boundary 


layer inside the nozzle! 


QUESTION: 


“How will SGS model affect our 
lipline RMS and farfield solutions” 
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Far-Field Comparison: PSD Spectrum at 100D from exit 
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Far-Field Comparison: Band-Limited OASPL (0.08 < St < 8.0) *8 
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SEM backup slide. See former presentation from last year 
when | added the SEM routines 


Add slide for different turbulent inflow methods. Comparison 
name of Spalarts SGT .... Maybe friction velocity plot. 
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EXACT LOCATION OF MICROPHONE POSITIONS 
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End cap averaging for FWH surface. Add more info in 
correspondence with Joseph. 
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